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Abstract. In the classical limit no work is needed to couple a system to a bath with sufficiently weak 
coupling strength (or with arbitrarily finite coupling strength for a linear system) at the same temperature. 
In the quantum domain this may be expected to change due to system-bath entanglement. Here we show 
analytically that the work needed to couple a single linear oscillator with finite strength to a bath cannot be 
less than the work obtainable from the oscillator when it decouples from the bath. Therefore, the quantum 
second law holds for an arbitrary temperature. This is a generalization of the previous results for zero 
temperature [HE]; in the high temperature limit we recover the classical behavior. 

PACS. 03.65.Ud Entanglement and quantum nonlocality - 05. 40. -a Fluctuation phenomena, random pro- 
cesses, noise, and Brownian motion - 05. 70. -a Thermodynamics 



1 Introduction 

The second law of thermodynamics [3] is considered one 
of the central laws of science, engineering, and technology. 
For over a century it has been assumed to be inviolable 
by the scientific community. Over the last 20 years, how- 
ever, its absolute status has come under increased scrutiny 
[4]. Challenges to the second law have recently attracted 
big interest with consideration of the miniaturization of 
a system under investigation, especially at low tempera- 
tures where quantum effects are important [4, 5j . In con- 
trast to common quantum statistical mechanics which is 
intrinsically based on a vanishingly small coupling be- 
tween system and bath ("thermodynamic limit") , the finite 
coupling strength between them in the quantum regime 
causes some subtleties that must be recognized. The quan- 
tum thermodynamic behaviors of small systems have the- 
oretically been investigated intensively and extensively [H 
and experimentally been examined [4,10, 



EI. 

The problem of a quantum linear oscillator coupled to 
an independent-oscillator model of a heat bath (quantum 
Brownian motion) has been extensively discussed [12,13, 
fT4jfT5|ll6| . The validity of the quantum second law has 
recently been questioned in this scheme at zero tempera- 
ture [5j fT7lll8j by the fact that the coupled oscillator has a 
higher average energy value than the free harmonic oscil- 
lator ground state, which would not be in compliance with 
the second law. By means of a cyclic coupling/decoupling 
process one might expect to extract useful work from a 
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single bath. However, this claim turned out incorrect; the 
apparent excess energy in the coupled oscillator cannot 
be used to extract useful work, neither for the well-known 
Drude damping model with a cut-off frequency for the 
spectral density of bath modes shown by Ford and O'Connell 
in pQ nor for both discrete bath modes and continuous 
bath modes with the generalized realistic damping mod- 
els [2], since the minimum value of the work to couple the 
free oscillator to a bath takes above and beyond this ex- 
cess energy. Therefore, the quantum second law for zero 
temperature is inviolate. 



In this paper, we would like to discuss the second law 
in the quantum Brownian motion at an arbitrary tempera- 
ture. We will obtain an analytic expression for the second- 
law inequality which can explicitly be shown to hold for 
the Drude model, which is the prototype for physically 
realistic damping. It is known [19] that a finite frequency 
cut-off reflects the physical fact that the bath cannot react 
instantaneously to a change of the system oscillator, and 
that in the absence of the cut-off, some observables such 
as the variance of the system momentum diverge. There- 
fore, the physically unrealistic cutoff-free damping models 
considered in [2] will not extensively be considered here 
(cf. Sect. [5]). From the result of this work, the appearance 
(or disappearance) of quantum effects versus thermal fluc- 
tuation over the different temperatures will also be seen 
explicitly. Let us begin with a brief review on the basics of 
the quantum Brownian motion. We will below adopt the 
notations used in |20"1. 
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2 Basics of quantum Brownian motion 

The quantum Brownian motion under consideration is de- 
scribed by the model Hamiltonian 



H = H s 



Hb-sb 



where 



V M 2 -2 

H s = — 1 uiq 2 

2M 2 ° q 



(1) 

(2) 
(3) 



The Hamiltonian H s -sb splits into the bath and the cou- 
pling term, 




N / -2 



2m, 1 2^ 



N N c 2 

-9 E c i £ j + ^ E - 



3=1 



- 2m,- u> 2 

j=l J 3 



(4) 



(5) 



From the hermiticity of Hamiltonian, the coupling con- 
stants Cj are real-valued. Without any loss of generality, 
we assume that 



Wl < U>2 < ■ ■ ■ < lun-i < wn 



(6) 



By means of the Heisenberg equation of motion, we can 
derive the quantum Langevin equation 

M 'q + M J* ds 7 (t - s) q(s) + M UJ 2 q = |(t) , (7) 

where we used p = Mq, and the damping kernel and the 
noise operator are respectively given by 



1 N c 2 



3=1 3 1 

N 



cos(wj t) ; £0) = -Af 7 (t) q(0) + 



^2 Cj l Xj{0) cos(Uj t) + h^L sin^ t) f- . (8) 

Here, (£(t)) Pb , = and (i(t) i(t')) Pb , = M 7 (t - f)//3 
for the initial bath state with the shifted canonical equi- 
librium distribution, py = e~^ Hb - ab ^ QS> /Z^ [19] where 

/3 = 1/ksT, and Zg b ' is the normalization constant (i.e., 
the partition function) . The Fourier-Laplace transform of 
lit) is [20] 



7(w) 



/W 



n 

E 



j j j 



(9) 



Introducing the spectral density of bath modes as a char- 
acteristic of the bath, 



N c 2 



3=1 



(10) 



3 ^3 



we can also express the damping kernel as 
2 f°° duj J{uj) 



lit) 



M 



■ cos(o;t) 



(11) 



We now consider a response function [20], xit) 
[q(t), q])a, where the expectation value (• 



with respect to the equilibrium state, pp 



is taken 



-0H 



with the partition function, Zp = Tre _/3ff . The Fourier- 
Laplace transform of xit) is then the dynamic susceptibil- 
ity 

XM = J?— 9 1 ■ -/ v (I 2 ) 

M UIq — LU Z — lLOj(iO) 

which plays important roles later. It is known [12] that 
the susceptibility x(w) can be rewritten as 



N 



IK- 2 



1 .7 = 1 



M ^ 



(13) 



k=0 



«2) 



where the normal-mode frequencies, Wfc of the total system 



if satisfy Wq — w 2 . — iujkli^k) 



generality, we here assume that 

Qq < uj\ < ■ ■ ■ < u)jv— i _■ 
It can then be found [Lip] that 



0. Without any loss of 
(14) 



(15) 

and wo < wo < wjv- 

The damping function 7 (w) in ([9]) can be rewritten as 



(16) 
(17) 



'V^y Af Jo 7T u,' ^o)M 
7(o> 



00 dui' J(^') 



P 



which resulted from the Fourier-Laplace transform of ifTT]) . 
Here, we used the well-known formula, l/(x + i0 + ) = 
Pil/x) — in5(x) for x = w' — w. Equation lfl7|) is con- 
venient for the case of a continuous distribution </(w) of 
bath modes; for the simple Ohmic case, Jo(w) = Mj q uj 
with an w-independent constant j , we easily have 70 it) = 
27 <$(£), and 7o(w) = 7 with a vanishing principal (or 
imaginary) part in fT7|) , while for the Drude model where 
Jdioj) = M j uj oj^/iu) 2 +w^) with a cut-off frequency u>d, 
we have 7 d (i) = 7 D w d e" 



7d(w) 



w 2 + w3 



and 



7o^d^ 
UJ 2 + UJ 2 , 



(18) 



Here, Jdito) behaves like Jo(w) for small frequencies (with 

LO d -> OO). 

The model Hamiltonian in ([1]) can also describe the 
classical Brownian motion (see, e.g., [21]). In considering 
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the quantum second law below, therefore its classical coun- 
terpart will be discussed briefly in comparison. 

One might think of using instead of H s b in ([2]) its ro- 
tating wave approximation H^' = K j (^J + ^ ty), 
which has for uo = u)j for all j, energy-conserving terms 



only. Here, a 



2h y ~ 



2h ""J ~ 

j== Pj, respectively. In this case, though, we cannot 

a / 27fi j Tito j 

observe any excess energy of the coupled oscillator at zero 
temperature since the system oscillator and all bath os- 
cillators remain unchanged in their ground states, respec- 
tively, with no entanglement (see also [22ll23j): from the 



Heisenberg equation, ihH. 
we can obtain 



(r) 



H 



H 



(r) 



i 
i 



[H s , H«] = [H 8 , H%>], 

Kj (a+Sj - a St) 

j 

''-Kj (a St - tils) , (19) 



which yield H^^i) = H^'\tpi) = for the initial state 

\ipi) = 1 0) |00 • ■ ■), respectively. For the full Hamiltonian H 
in (!]), on the other hand, we have, after a fairly lengthy 
calculation, 



rWi 



(a - at) (Sj + St) + 



2iM 



t^2 



^21 



/ ^ 2 ' 

7 J J 



(20) 



which clearly gives rise to 7^ 0. 



3 Formulation of the quantum second law 

From the fluctuation-dissipation theorem [19] , we can eas- 
ily have 

\{q{t x )q(t 2 ) + q(t 2 )q(ti))p = (21) 
I / °° du, coth (^) cos{a^ 2 - h)} lm{x(co + 1 0+)} , 

which immediately yields 

(? 2 >/3 = I Jo°° <*" coth (^r) Im {x(^ + * 0+)} (22) 

= £ / °° dw c 2 coth (^) Im{x(o; + 1 0+)} (23) 

and thus the energy of the coupled oscillator 

E S (T) := <F S ), = /"^(yg + ca) coth 

xIm{x(w + «0+)}. (24) 

In comparison, the internal energy of an uncoupled (or 
free) oscillator is [3] 

e(co,T) = ^0(5 + = ^ coth ^' ( 25 ) 



where the average quantum number (n)p = l/(e l3hu "' — 1 
Its classical counterpart appears as e c ; (T) = 4 . With ([25 
equation (|24|) can now be transformed to an expression 



E S {T) = E s (0) + AE S (T) 



(26) 



where 



Mh 



du> 



AE S (T) 



Mh 

2jri 



du (lu 2 + uj 2 ) x(u) (n)p • (27) 



Here, the integration path is a loop around the positive 
real axis in the complex w-plane, consisting of the two 
branches, (00 + ie, ie) and (— ie, 00 — ie) [24]. Therefore, 
E S (T) for the discrete bath modes can be exactly obtained 
in closed form from the residues evaluated at all poles {wj,} 
of x(oS) in lfT3|) on the positive real axis. Then, equation 
(l26l reduces to 



1 N { 
E s(T) = -]>>(^,T) 1 

fc=0 I 




To study the quantum second law below, we need two 
different Helmholtz free energies for the oscillator coupled 
to a bath. 

We first consider the Helmholtz free energy of the cou- 
pled oscillator, F S (T) = E S (T) — TS S [3] with its entropy 
S s = —dF s /dT. This can easily be solved for F s such that 



pT E 3 (T') 
To 

By requiring that the entropy 

S s (p) = k B 0E s (p) 



F S (T) = Tl-ti^dT 







C 



k B / E s (P')dp' 

00 



(29) 



C (30) 



with /3q = l/fceTo vanish at zero temperature, we can 
determine the constant C of integration; using ([24]) we 
easily obtain 



C/k B 
where 



Mh 
2 



£ / °° dcj (c 2 + a, 2 ) Im{ X (u; + l 0+)} A{u>) , 

(31) 



coth 



L„ 43' coth 



13'huj 
2 



(32) 



The asymptotic series, cothz = 1 + 2^^ 1 e 2lz then al- 
lows equation lj32j) to become 



A(u) = f° df3' coth ^ 



(33) 



which clearly makes equations f29|) and |30)l , respectively, 
independent of To (or and accordingly uniquely deter- 
mined. 
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The other Helmholtz free energy T s (T) needed for the 
second law is described as follows; the minimum work 
required to couple a system oscillator at temperature T 
to a bath at the same temperature is equivalent to the 
Helmholtz free energy of the coupled total system minus 
the free energy of the uncoupled bath [12] . This minimum 
work can then be obtained as the free energy T S {T) = 
— A InZp, where the canonical partition function Zp = 

Tr e -1311 /Tib e~^ Hb . Here, Tr& denotes the partial trace 
for the bath alone (in the absence of a coupling between 
system and bath, this would exactly correspond to the 
partition function of the system only). By means of the 
normal-mode frequencies Qk, we easily get 



A? 



ft e -/3tofc(n fc +^) 



z fJ = 



k=0 n k =0 



N 



(34) 



[ e -^i n i+i) 

j=l nj =0 



which yields 



T S (T) = ]T f(^,T) - ]T /(^,T) (35) 



fe=0 



with the free energy of an uncoupled oscillator 

hjj l 



f(",T) 







In (1 - e- im ") 



In f3htu 



(36) 



with 



The classical counterpart of l|36p is / c i(w, T) 
h <C 1. Equation (|35| can then be rewritten as [12] 

1 f°° id 
T S (T) = - / du>f(w,T)1m< — lnx(w + iO+) 
tt J [duj 

(37) 

in terms of the susceptibility. We note here that for an 
uncoupled oscillator, Im { ^ In \(ui + i0 + )} — > n 5{uj — 
ojq) and thus T s (T) -> f(uj ,T). Similarly to (37|), we 
can also obtain the energy required to couple a system 
oscillator to a bath, 



£.(T) = £e(w fc) T) -5>K,T) 



(38) 



fc=0 



= -/ duje(uj,T) Imj-f- lnx(w + i0+) j . 

7T J rfw J 

From e(w,T) > f(u),T) (the equal sign holds for T = 
only) with the frequency relationship in (fl5| , we can easily 
get £ S (T) > T S {T). 

We now consider a cyclic process composed of the cou- 
pling of a harmonic oscillator to a bath and then the 
decoupling of the oscillator from the bath (the coupling 
constants Cj — ► 0). The free energy change on comple- 
tion of the coupling process is J- S (T) — f(u>o, T), whereas 
the maximum useful work obtainable from the oscillator 
only in the decoupling process is the free energy difference 
F S (T) — /(woi T) which cannot be greater than the energy 



change E s (T) — e(u>o, T). Here it is assumed obviously that 
the extraction of energy from the bath is impossible. The 
second law can then be expressed as an inequality 

T S {T) - f(uj ,T) > E S (T) - e{io a ,T). (39) 

In obtaining lj39j) we used the conceptional difference be- 
tween J- S {T) and F S (T) ("operational asymmetry") |26j. 
For zero temperature, this inequality, obviously, reduces 
to T S {Q) > E s (0), the validity of which has been explic- 
itly proven for the Drude damping model pQ and for the 
discrete bath modes, by means of E s (0) in ((28)1 and ^ s (0) 
in ([35)1 with the frequency relationship (fl5| , and the gener- 
alized realistic damping models of continuous bath modes 
[2]. For non-zero temperatures, on the other hand, it is 
very non-trivial to investigate the validity of inequality 
(f39"l) with E S (T) and J- S (T) for the discrete bath modes. 
For the continuous bath modes, the evaluation of E S {T) 
and T s (T) clearly depends on the parameters of the damp- 
ing model considered. We will below discuss inequality 
(j39"l) explicitly within the Drude model which is the pro- 
totype for physically realistic damping. 



4 Discussion of the second law within the 
Drude model 

It is convenient in the Drude model to adopt, in place of 
(u)q, uid, jo), the parameters (wo,/2, 7) through the rela- 
tions p] 



^0 : = w o ITT^ ' Wd := Q + T 

_ f2(J3 + 7 ) + w^ 
1° •— 7 (f2 + 7 ) 2 



(40) 



Substituting equation 1|T8|) with lj40j) into {T2JI , we obtain 
the susceptibility 

~ , \ 1 u + i{f2 + zi + z 2 ) 

M (uj + ill){u) + izi)(u> + 1Z2) 

where Zi = j/2 + iwi and z 2 = j/2 — iwi with wi = 
a/\Vq — (j/2) 2 . First, we consider the overdamped case 
(7/2 > wo), where Zi, z 2 > 0. We then have 

z=i — 
where Wi = f2, UJ2 = %i, u>3 = z%, and the coefficients 



A 



(i) 



Zl + Z2 . \ 



(2) 



(r2-zi)(z 2 - ft) ' 



r? + z 2 

(zi - ft)(z 2 - zi) 



A 



(•■i) 



j? + zi 



d " (z 2 - fi)(zi - z 2 ) 

Here, we note that 

3 3 



(43) 



£A« = ; ^ ^ 2 = 0. (44) 
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To obtain an explicit expression for the energy Eg (T) 
of the coupled oscillator, we first consider the integral in 
f2Tj) : by performing a contour integration with the aid of 
|42|) and the identity 



coth 



(59 -£('♦'£* 



to 



where v n = 2irn/f3H, we can have 

3 ( oc 
~JM E X d { ^TJT + 2 E ~ v* 



, (45) 



(46) 



Z=l 



n=l 



With ft — > 0, this reduces to its classical counterpart, 



relation [20] 



(q(0)q(t) + q(t)q(0))p = -**<$(0)$(t) + g(t) q(0))p 

(47) 

it can eventually be found that 



(0 



(48) 



E 



With ft — > 0, its classical counterpart easily appears as 

\ = e cl (uj ,T). We can 



p(d)frp\ _ 1 \ (Q _l -^Q _ 1 



also express (|46|) and ([47)) at t = 0, respectively, in terms 
of the Digamma function |27| 



dy 



OO - OO 



-c e + 



— (49) 

n=l n=0 y 



with the Euler constant c e as 

3 



(50) 



i=i 



2=1 

Here, we used (|44|l . Equation (|48|l can thus be rewritten 



^(T) = |E^(-o 2 -^ 2 ){^ + ^(^)} 



i=i 



(52) 

With the aid of the asymptotic expression, ip{y) =lny — 
~ky ~ S^Li 2nP" wu ^ n Bernoulli number B n [27], the 

zero-temperature value £^'(0), clearly, reduces to (|66p 
derived in [2] (see Appendix). 



Let us now consider the free energy Fs(T). By sub- 
stituting {41]) into {37J with the identity, ln(l + y) = 
~ Y^=i {~y) n / n > we can easily obtain 



1 00 1 

Ti d \T) = J*0(0) + — ^ - A»C8) , (53) 

n — 1 



where 



wi In 



7/2 - Wi 
7/2 + Wi 



(54) 



with wi = yj (7/2) 2 — Wq, as was derived in [2], and 



4.(0 = E W / 

.._n JO 



, x r°° c—nPtiUnV 

Cm) ' ' e 

3 



dy- 

2T + 1 



(55) 



^ {sin(n/3ftov) Ci(n/3ft^) 
(n(3tkOy) si(n/3fta^)| 



£t=0 

cos 



(56) 



with = 1, t^ 1,1 = = = -1, and ^0 = u d (see 
also Appendix). By the substituting the classical quantity 
fci{uj,T) into |37|) with J* dy/(y 2 + 1) = arctanx and 
J °° dyhxy/(y 2 + 1) = [28], we can also obtain the clas- 
sical counterpart c iTs d \T) = —jg^-^ + fcl(^o,T) = 
f c l(ui Q ,T). 

For the underdamped case (7/2 < wo), equations (|52"ll 
and lf53j) can be found to hold as well, respectively, being 
expressed in terms of the functions with complex-valued 
arguments. By showing the validity of inequality (|39| for 
underdamped and overdamped cases, 

K d {T) := JFW (T) - /(wo, T) - E^(T) + e(u; ,T) > 

(57) 

as in Fig. [H we see that there is no violation of the quan- 
tum second law; in fact, Kd(T) vanishes asymptotically 
with the increase of T. 

Comments deserve here. In the classical treatment both 
sides of inequality i|39j) vanish, namely, c iKd(T) = 0, which 
clearly means that no work is required to couple a linear 
system to a bath at the same temperature, and no energy 
change in the system is obtained during the decoupling 
(and the coupling). In the quantum treatment, on the 

other hand, Ei d \T) and J-s d \T) depend on the damp- 
ing parameters, respectively. Therefore, while both sides 
of ([39| become vanishing in the high temperature limit 
(equivalently, ft — > 0), they actually do not vanish espe- 
cially in the low temperature regime. This non-vanishing 
behavior stems from the system-bath entanglement in- 
duced by the finite coupling strength between them, which 
leads to the deviation from = e~ t3 ^ B /Zi 8 ' for the 

reduced density matrix p s d '(T) being, clearly, damping- 
parameter dependent. 



.(1) _ J2) _ (3) _ 
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In fact, we have Kd{T) > especially in the low tem- 
perature regime (see Fig. [J). This strict irreversibility over 
a single cycle composed of the coupling and decoupling 
process appears from the fact that the system-bath en- 
tanglement induces the entanglement between any pair 
of infinitely many bath oscillators ("entanglement swap- 
ping" [29]), which cannot completely removed over the 
system-bath decoupling process. Therefore, we essentially 
cannot recover the original state of the bath, pi ^ and thus 

that of the system, p^K As a result, (T) — f(aj ,T), 
being the minimum work required for the entangling in 
the coupling process, is greater than the energy change 
Ei d \T) — e(ujQ, T), which can necessarily not be less than 

the free energy change, Fg d \T) — f(uj , T) being the max- 
imum useful work obtainable from the system only in the 
decoupling process. With the increase of T, however, the 
strict irreversibility shrinks (Kd{T) — > + ) since the ther- 
mal effect dominates the quantum effect. In the classical 
case, on the other hand, this operational asymmetry, intro- 
duced in the last paragraph of Sect. [3] disappears at an 
arbitrary temperature, namely c iTs d \T) = c iF s (T) = 
fa(u ,T). 



5 Comparison with the Ohmic model 

Let us briefly consider the Ohmic model (as a cutoff-free 
damping model) for an arbitrary temperature to compare 
with the Drude model considered in Sect. [4] For zero tem- 
perature it is known that K o (0) = JT\o) - E [ s o) (0) 
vanishes, where 



(o) 



diverges logarithmically, while its Drude-model counter- 
part Kd{0) — > Egj/nwo in the limit ujd — * oo (equiva- 
lent^, i? — ► oo) where E g is the ground state energy of a 
free oscillator. 

For the overdamped case (7 D /2 > u>q), the susceptibil- 
ity in (fl2|) appears as 



1 



(59) 



M (u> + iuji) {u> + iu)2) ' 

where u)\ = j /2 — w > and 1V2 = 7o/2 + w > with 
w = ^(7 /2) 2 - u%. Substituting ((59]) into l[2T]). we can 
obtain 



(60) 



j=l { n=l J 

with A^ = — 1 and Ao = 1 [30], from which, similarly to 
El, 



~2\(o) 



~ 2 



Substituting {60]) into (|47|) . we can also get 



yi h ~ 2w 



3=1 



^ ^ n 



(62) 



which diverges logarithmically for an arbitrary tempera- 
ture (note that linin->ooQ^fc=i k ~ mn ) = 



compare 



this with (q 2 )^ in ([5l]l being convergent. From ([BTjl and 
{62]) , the energy of the coupled oscillator is 



e { s °Ht) = & £ A« {(«g - ^) + |v, } 



which clearly diverges. With 7l — > 0, its classical counter- 
part, however, reduces to c iEs°\t) = e c i(uj ,T) (it can 
also be found [9] that Es(T) is identical to the energy 

£s°\t) obtainable from {381) )• Further, similarly to {53]), 
we can easily obtain 



T { °\T) = Ti°\0) + AF s (o) (T), 



(64) 



where 



11 

ZiJTW(T) = - ^ £ - £ { shiMH') Ci(n^) 
^ n=l " 3=1 
— cos(n(3ftwj) si(n/3hujj)} . (65) 

The free energy Ts°\t) also diverges. With ft — > 0, we 

get c iT ( s o) {T) = f c i(uj Q ,T). Equations {51 and {64]) can 
be found to hold, respectively, for the underdamped case 
as well. 

For comparison with the Drude model, we take the 



limit Q {or u>d) — > oo in (|43| so that A 



(i) 



0, A 



(2) 



(3) 

1/(2! - z 2 ), and A^ J; -> -AJf\ Then, it can easily be 

shown that E<f\T) v> £i o) (T) and ^ s (d) (T) -y» ^ s (o) (T); 
the second term on the right hand side of {53]) reduces 
to AF s (o) (T) in {64]), however, the first term ^ d) (0) -n 

Ts°\o) as discussed in [2]. As a result, the second-law 
inequality (|3"9j) for the Drude model with u>d — > oo is not 
equivalent to that for the Ohmic model. Whereas the clas- 
sical counterpart c iK (T) = 0, both sides of {39]) come to 
diverge differently so that it is non-trivial to explicitly 
evaluate K a {T) for this unrealistic damping model. 



l(2) 



6 Conclusions 

In summary, we have studied the second law in the scheme 
of quantum Brownian motion at an arbitrary temperature. 
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It is clearly a generalization of the previous works for zero 
temperature by Ford and O'Connell [I] and by the authors 
of the present paper [2]. It has been shown for the phys- 
ically realistic damping model that the work needed to 
couple a system oscillator to a bath at the same tempera- 
ture cannot be less than the work obtainable from the os- 
cillator only when it is extracted from the bath; especially 
in the low temperature regime the apparent irreversibil- 
ity, Kd(T) > 0, stemming from the system-bath entangle- 
ment was found, which is different from the behavior of its 
classical counterpart, c iKd{T) = 0. Therefore, the quan- 
tum second law holds for an arbitrary temperature. The 
question about the validity of the quantum second law for 
a broader class of quantum systems than the quantum 
Brownian motion considered here, especially non-linear 
systems coupled to a bath, clearly remains open. 
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A : Mathematical supplements 

It has been shown in [2] that for the overdamped case 
(w < 7/2), 



E { s d) (Q) 
where 



2tt 



{A(w ,/2,7) + S(w ,r2, 7 )} , (66) 



(wg+r? 2 ) in 7 2 /4 - rtwg - 7 /2) + 7 a / 



(67) 



^7/2 + wi J 



with wi = \J (j/2) 2 — Wq , and 



*(wo,fl,7) = ( ^7Hw|- 2 ^:^) Mfl/wo). (68) 



In derivation of equation (|56|) from (|55| , we used [28j 
j^dy^^j = \ {sin(a6) Ci(a&) - cos(a6) si(afr)} , (69) 

where a, b > 0; the sine integral si(y) = — dz sm W — 
-f + Si(y) with Si(y) = tfdz 



y a-, sin O) 



tegral Ci(y) = - dz 
with the Euler constant c P 



o 



and the cosine in- 



= c e + In y + J y dz 
0.5772156649- • •. 



cos(z) — 1 
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Fig. 1. 



Fig. HJ y = Kd{T)/hwo versus temperature T; from bot- 
tom to top: (J? = 1 and 7 = 3/2 underdamped) , (J? = 1 
and 7 = 4 overdamped), (42 = 5 and 7 = 3/2 under- 
damped), and (f2 = 5 and 7 = 4 overdamped); here, 
h = ks = wo = 1; in the high temperature limit, we 
have the classical behavior, y — > + . 



